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On certain Possible Cases of Steady Motion in a Viscous 

Fluid. 

By Thomas Craig, 

Johns Hopkins University and United, States Coast and Geodetic Survey. 



The following paper contains, first, some general principles governing steady 
motion in viscous fluids ; second, the detailed working out of two cases, (i) a 
sphere moving with constant velocity in the direction of the axis of x, (n) an 
ellipsoid moving uniformly in the same direction. The results obtained are 
certain to hold for slow motions, though they have been obtained without that 
assumption, but it is not proved that the prescribed conditions will exist for rapid 
motions. If it can be shown that a velocity can be chosen for the moving body, 

so that the quantity 

A 2 u . dx -\- A 2 v . dy + A 2 w . dz 

shall be an exact differential, then the solution below given will hold for that 
case, and for that case only. 

Part of what immediately follows I have already given in another place, but 
it is repeated here for convenience. 

The expressions for the fluid pressure in different cases are given in the 
" Journal of the Franklin Institute," for October, 1880. The values found for 
the velocities of a fluid particle when a sphere moves in any direction in the 
fluid are given in the " Philosophical Magazine " for November, 1880. A slight 
error exists in these values as there given, which is corrected here. 

Denote by u, v, w the component velocities of a fluid particle in the direction 
of the axes x, y, z ; p, the density of the fluid at the point x,y,z; p, the constant 
density, and [i, the coefficient of viscosity ; the kinematic coefficient of viscosity 

or the ratio of p. to p will be denoted by k, i. e. - = h. 

The equations of motion of an incompressible viscous fluid are now 
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du „. 1 dp . . .„ 



%=r-\% + ^\ CD 



dw __ Idp . 



X, T, Z being external forces. In what follows we will suppose the forces 
X, Y, Z to possess a potential. If the motion of the fluid is caused by a body 
which has been projected in it, and is acted upon by forces due to a potential, 
the potential must be of the form 

Ax + By + Cz, 

for the motion to be steady relatively to the body, as a constant resistance has 
then to be overcome. If the body is at rest and the liquid streaming past it, 
the potential must contain a term of the form 

Ax + By + Cz 

at infinity, to keep up the steady motion ; otherwise the motion would die away 
and the liquid come to rest from the presence of factors of the form e~ pt * 

Denoting by £ 77, £ the component angular velocities of the fluid particle at 
the point x, y, z, we have 

j. 1 /dw dv\ 

±~2 \dy ~ dl)' 

_ 1 (du dw\ 
V-2\^~^)* ( 2 ) 

1 /dv du\ 

±~2 \dx ~ ch}J ' 

From these follow readily the known relations 

^ = »(§-£> (s) 

*»- »(S-2> 

* I am indebted to Mr. Greenhill of Emanuel College, Cambridge, for the above remarks, and also for 
many other most valuable suggestions. 
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Denote by D, the resultant angular velocity ; then 

2 = ^ + 7? a +£ 2 ; (4) 

also write 

2 q = u 2 + v 2 + w\ (5) 

The internal friction involves a certain dissipation of energy ; the function 
expressing the rate of dissipation per unit volume has been called by Lord Eay- 
leigh the " dissipation-function ; " denoting this by E, we have 

_ rt ( /du\* , /dv\ 2 . /dw\2 . 1 /dw , dv\ 2 

(6) 

. 1 /du dw\ 2 1 /dv . du\ 2 ~i 
~1~2 \Jte + fa) + 2 \dc~T~fy) y 

This can be given in a different form by obtaining the expression for A 2 q ; this 
is readily found to be 

A 2 q — uA 2 u + vA 2 v + ioA 2 w 

du\ 2 f^y, /*A 2 

dx) "+" [dyj ' [fej 



+ (£)*+ (i) ■+ (sy (7) 



fdw\ 2 /^Vi /^w\ 2 
+ ^J + ^J + ^&j" 



^ 



Add 2fi 2 to g - , and we will eliminate terms of the form 



dw dv 

dy "dz-> etc ' 

Then compare the resulting form of equation (6) with equation (7), and we 
have at once 

E = 2^{A 2 q — (uA% + vA 2 v + wA 2 w) - 2 O 2 } (8) 

or 

uA 2 u + vA\ + wA 2 «o =A 2 q — 7^E—2Q?. (9) 

In equations (1) the quantities on the left-hand sides may be replaced by 

du . du . dv . dw . . 

Tt+ U fa + V dx + W fa+ 2 ( W V- v 0> etc., 
or by 

Tt + £ + 2 ( Wr ? ~ w 0> etc. 
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On making these changes, the equations of motion become 

du , d(V+q) 1 dp n _, A . 

* + ^P +p I + 2 «"^) = *^^ (10) 



P 



¥ + -V-+;s + 2 («f-«i) =3 ^ 



To these is to be added the equation of continuity, 

i+|+i=o. (ii) 

Write 

then, introducing the conditions for steady motion, (10) become 

^ + 2(w v -vQ = kA\ 

%+2(ut-v>£) = kA% (12) 

dP 

■& + 2 K - uy) = *A 2 «>. 

If we assume that the quantities A 2 u, A 2 v, A 2 w are the first differential co- 
efficients with respect to x, y, z of a function $, these equations become 

d(P-kQ) _ 

Tx = -2(w v -v£), 

d(P-kQ) _ 

Ty = — 2(uC-wg), (13) 

d{P-kQ) _ 

Jz = -2(vi-urj). 

Multiplying these by u, v, w respectively, and then by £, rj, £, and in each case 
adding the results, we have, writing for brevity 



@ = P — JcQ, 

= 0, 

(14) 



d® , d® . d® 
U Tx + V d-y+ W Tz = > 



£ d® d® d® _ 

$ dx"T~ V ty ^~ £ H ~ U ' 
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and also 

d® ...» ._ „. 

^ = 2*2 sin S, (15) 

where q = V2~q is the current velocity, and 8 is the angle between the stream 
line and the vortex line at the point x, y, s. Hence the conditions that the state 
of motion of the fluid for which 

A 2 u . dx + A 2 t> . dy + A 2 w . dz 

is an exact differential, are as follows : It must be possible to draw in the fluid a 
system of surfaces, © = const., infinite in number, and each of which is covered 
by a network of stream lines and vortex lines. This is the property denoted 
by equations (14). The product q'D, sin Bdn must be constant over each such sur- 
face, dn denoting the length of the normal drawn to the consecutive surface of 
the system. These results are identical in form with those given for a perfect 
fluid by Professor Lamb in his work on Fluid Motion. In order that 

A 2 « . dx + A 2 v . dy + A 2 w . dz 
shall be an exact differential, the equations of condition 

A 2 £=0, A 2 t ? = 0, A 2 £=0, (16) 

must hold. 

If Ave assume that the motion of the fluid is so slow that squares and pro- 
ducts of the velocities may be neglected, equations (1) become, when a potential 

exists, 

M 2 * = f> (17) 

7A2 dU 

M 2 W = w > 

where 

u= r+ff> (18) 

In this case the quantity 

Ah( .dx.-\- Ah . dy + A 2 w . dz 

is obviously an exact differential, and the equations of condition 

A 2 £ = 0, A 2 V = 0, A 2 £ = 0, 
are satisfied. 

From equations (12) we have, in every case, 
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« (g - MS,) + v (f - MS.) + to (f - MS„) = 0, 

f (f - MS,) + , (f - MY.) + I (f - MSo) = 0, <19) 

so that the conditions for steady motion hold now as in the particular case just 

mentioned ; but in this case the surfaces in the fluid are given by the differential 

equation 

dP — Tc (A 2 u . dx + A 2 v . dy + A 2 w . dz) = 0. (20) 

Write for convenience 

L = vC — wrj, 

M=w£-u£, (21) 

iV= Urj — v£. 



Equations (12) now become 



~-2L=kA\ 

~-2M=kA\ (22) 

^-2JST=kA 2 w; 
and from these, by differentiating for x, y, z, respectively, and adding, we have 

A'i> = 2(g + f + f). (23) 

The same equation holds when the motion is not steady ; for if we differen- 
tiate equations (10) for x, y,-z, respectively, and add, the terms containing the 
differential coefficients of u, v, w with respect to t will disappear by virtue of 
the equation of continuity. Integrating (23), and substituting for P its value, 
there results 

? = Hn///(£+f +£)*** + "+.}■ (24 > 

In the general case G is a function of the time ; for steady motion, however, 
it is a constant. The quantities L, M, N may be the first differential coefficients 
with respect to x, y, z of function of x, y, z ; suppose such a function W to exist 
that we have 

2M=%> (25) 
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equations (22) become in this case 



d(P-<a) 

fiT 2 = ^ V 



d -^ = M\ (26) 



dy 
d(P-&) 



and from these- results 
and consequently 



A 2 (P-^) = 0, (27) 



in which 

A l = dx'* + ^" 2 + <fe' 2 ' ( 29 ) 



Equation (24) can be thrown into another form by very simple transformations. 

Write 

Z> 2 = £ 2 + IP + JSP; (30) 

then, denoting by a, f$, y the direction-cosines of the vector D, 

X= aZ>, 

M=0D, (31) 

N=yD. 

Substituting in (30) the values of L, M, JSF, we have 



^=^[i-(?l+^+S) ! } < 82 ) 



or 



Z> = g'n sin 8, (33) 

where 8 is the angle between the stream line and the vortex line at the point 
x, y, z. Now let a, b, c denote the direction-cosines of a normal to the closed 
surface containing the fluid ; then 

aa + b/3 + cy 

is the sine of the angle between the normal and the plane containing the instan- 
taneous axis of rotation O, and the direction of the resultant velocity q, or 

sin <f> = aa + bfi + cy. (34) 

Similarly, if we denote by a', b', c' the direction-cosines of the line joining (x, y, z) 
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to (x', y', z') and <f> the angle between this line and the above-mentioned plane, 
we have 

sin f = a'a + &'/8 + c'y. (35) 

Take now the triple integral in (24), 

f dL dM dN\ dx'dy'dz' m 



1 rrr/dL dM dN\ 
2* JJJ \dx> + dy + dz'j 



2 
this is 

= 2? J J ( ai + 6Jf + ciVr ) "7 + £ JJJ ? «fe%W, 

where Jo- is an element of the bounding surface. By virtue of the above equa- 
tions, we have 

aL + bM + cN = (aa + bfi + cy) D = gfil sin 8 sin <£, 

and similarly 

a'L -j- b'M + c'N = g'ft sin 8 sin </>' ; 
therefore 

'rfi rfJIf dN\ dx'dy'dz 1 



2-n 



(36) 



1 CCCf dL — , d:N \ dx ' d i 

1* JJJ \dx> + dy 1 + dz'J r 

1 P Pq'U sin 8 sin <f> , , 1 f /* /* o'O sin 8 sin <£' ,,,,,, 
= 2^JJ r ** + &]]] ? dxdydz, 

and finally (24) becomes 

| = # _ (F+ 2) + _ J J do- + ^JJJ ? dx'dy'dz'. (37) 

If the motion in the fluid is a screw motion, i. e. if the direction of motion be 
along the instantaneous axis of rotation, we shall have 8 = 0, and conse- 
quently 

l p = G-(V+q). (38) 

If the plane containing the direction of motion and the instantaneous axis of 
rotation be always normal to the bounding surface, we shall have <f> = 0, and 
then 

t = a - ( v+ g) + 1 /// *» "°j ""*' uau. (39) 

We will pass now to the consideration of one or two particular cases. As- 
sume, first, that u, v, w are given by the equations 
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u 




dW 
dy 


dV 
~ dst' 


V 


= 


dU 
dz 


dW 

dx ' 


w 




dV 

dx 


dU 
dy' 



(40) 



The functions U, V, W must, as is well known, satisfy the equations of con- 
dition 

A 2 £7=-2£, A 2 F=-2t ? , A 2 TF=-2£, 

<^,<WdW_ ( 41 ) 

dx + dy + dz ~ °" 

Instead of the three functions U, V, W, we may introduce a single function $ 
and write 

TIT <£* d$ 

w= y dx-- x dy-> 

these quantities will satisfy equations (41), and give us for the values of g, tj, £, 

Z=2\ Z dy—ydz-)> 

1 / dtft d$\ 

Z = 2\yT X - X dy)> 

where 

<j> = — A 2 3>. 

The function <f> must also satisfy the equation 

A 2 <£ = 0, (45) 

since 

A 2 £=A 2 tj = A 2 £=0. 

For (f> we can take any homogeneous function of the n th degree satisfying (45) 
and then will have 

*§+?f+*t=»*> (46) 
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or (f> represents a solid spherical harmonic of the degree n. For the values of 
u, v, w we have now 

d l , <7* c?* , d<& ) , 

u = Tx\* + x te + y^ + z di;\+ x <t>> 
w = & J* + * & + y % + 2 & | + ^ 

Reverting to equations (3) for a convenient form of obtaining the values of A 2 **, 
A 2 «, A 2 w, we readily find for these quantities the values 



a2 d $ 



A 2 * = n f y , (48) 



4,, ~= ft 
a2 _ ^ 

A 2 w = n Tz 



The function $ of equations (13) is now the function n<f>, or Q is in this case 
a solid spherical harmonic of the degree n. Equations (40) involve the assumption 
that the motion is purely of a rotational character. If we for a moment abstract 
the friction in the fluid from consideration, the motion, if caused by a solid mov- 
ing in the fluid, will be irrotational, and therefore subject to a velocity potential, 
say \fj, satisfying the equation A 2 »/f = 0. We can then in general write the values 

of u, v, w in the forms 

_ <fy dW dV 
u ~ <fo "■" dy ~ dz' 

<U , dU dW 

v = -dy + -dz~-W (49) 

__ dp dV dU 
w ~ dz + dx dy' 

Equations (43), giving the values of the rotation components, will be, of course, 
unaltered by this change in the values of u, v, w, and we shall have, instead of (47), 

d ( , . . cfo . d<l> , d®) , 
U = crx\-* + ® + X dx- + ydy + Z dz\+ Xt t>> 

W -dz\-t + ® + X ^ + ydy- + *az\+Z<t>- 
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From equations (43) we derive at once 

^ <fy <fy 

^ + V^ + ^Tz~ ' ( 51 ) 

from which it follows that the vortex lines lie on the surfaces given by the equa- 
tion (j> = const. From equations (14) we have, however, since in this case 

A 2 u . dx + A 2 v . dy + A 2 w . dz 

is an exact differential, 

t d® d® , d® __ 

Sfa+Vdi/ + {>dz — U > 

and consequently the vortex lines in the fluid lie at the intersection of the sur- 
faces 

d> — const. 

(52) 

® = const. 

The surfaces © are fixed in the fluid, but the surfaces <p may move; their 
motion, however, will always be in such a manner that the above condition shall 
be satisfied. Equations (50) can be thrown into a simpler form by the following 

considerations. Write 

r = (a? + y 2 + z 2 ) h 



and 
Then 


log r — \. 

d$ d$ d$ d<b 

x d^ + ydy +z di ==r d7' 


but 


d Id 




dr r dk' 


and consequently 


<?* d& 
r efr ~dk' 


We have then 


u = ?dk(®- x l' + dT) + x <f>> 




y . d / . . d$\ . 

v = s>dk{- x i> + * + dk) + y ( i > > 




z d ( d®\ 



(53) 



The solution of the problem when the motion of the fluid is caused by a sphere 
moving through it is quite simple* We have first to determine the velocity 

1 ■' — - c : — = - — — : ' — ' ■ — - — ■ — — ; — =^ K — 

* See an article on this subject by the author in the " Philosophical Magazine " for November, 1880. 
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potential. In Lamb's " Treatise on Fluid Motion " he gives a solution, due to 
Stokes, of the problem of a sphere moving with uniform velocity in a viscous 
fluid in the case when the motion of the solid is along the axis of x. and the 
motion of the fluid is symmetrical around this axis. Special polar- co-ordinates 
are employed in obtaining the required solution, but from the general values 
above given for u, v, w, we can readily obtain the same results in a very simple 
manner. We will consider this case for a moment, as the forms of u, v, w thus 
obtained are of use in another and rather more difficult problem. At an in- 
finitely great distance from the origin the fluid is streaming along the axis of x 

with a velocity = — X, or 

u = — X, 

v = 0, 

w = 0, 

for all of the motion which is due to a velocity potential. We have then, as 
indeed we know from other considerations, 

$ = — Xx. 
Our former value of <j> was 

CO 



where 

this will now reduce to 

<f> = <j>i = — Ax. 

For the general value of u we take into account the friction terms. The quan- 
tities Li and i? 4 will all disappear with the exception of L x and R x , and these 
also vanish at infinity. For these we have the values 

a— i +■!©•. 

For u we have, then, 

which, on substitution of the values of L x , R x , fa, becomes 

[V 3 a 1 /a\ 3 "l 3X fa a«"| . 



or u 

But 

and 

Writing then 
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— \ fl _ ? -1 _ ?^? ^ _ ^ f ^ - 3 as 3 ~| 

— X L 1 I rj 4 r 3 4 |_ r 5 J* 



we have finally 



or 



Similarly, 



r'-Sx* __ d? 1 
r? dx 2 r 



xl_ d_ 1. 
I s ~ x dx r' 



S\a 1 v 



Xa S 1 y 



u — X — Xi + %fa + a^> 



v= l( a;Xi+ £)' 



to 



= I (^ 



+ 



^?Y 

dx J 



The determination of the resistance experienced by the sphere, supposing it 
to move along x with velocity = + X, is the same thing as the determination of the 
pressure upon the sphere supposed at rest and the fluid streaming past it with 
velocity = — X. The latter case is the one that we are dealing with, and to 
obtain this pressure we use the dissipation-function, employing the form 

E = 2 p {A 2 ? — (uA 2 u + vA 2 v + wtfw) — 2 XI 2 }. 
We had 

2 ? = tt 2 + ^ + «) 2 , 
then 

2A 2 g = A 2 (« 2 + « 2 + w 2 ); 

substituting the values of u, v, w, gives 

2 , = (X - 2 V + 2 (X - 2 X.) f + (£) + (|) + (f )", 

where for brevity we have written 
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to this add 

P — 4 ' 

now, introducing the values of Q and X lf we find readily 

~ a„ ,8^ -12 a!* , „„ 36r 2 + 72a; 2 , 4r 2 -60a: 2 , . S^-r 2 
2A 2 ? = a 2 ? + 0* J + a/3 ? + 4 aX — ^— • 



We also find easily 
and 



«/ ao , ao ■ ao x ,4^ — 20a; 2 , rt 12a: 2 + 4r> . . „ 3a; 2 -r 2 
2 (wA 2 « + «A 2 + wA 2 w) = a 2 ? h a/3 p \- 4aX — — 



4f2 2 = 4a 2 
Combining all of these, we obtain 



V- j ,12a; 2 36r 2 +72a; 2 36a: 2 ) 

— ^ A* ] r* ' r 8 ^[r 2 A 4 ~ j 35 •" A* r 10 



=i2 M j5-¥i , s;+ / .^}i-si. 

Writing 

- = cos 5, 

multiplying E by 2 wr 2 sin J0dr, and integrating from 6 = to 9 = ir and from 
r = a to r = ao, we have, for the total rate of dissipation of energy, 



(3/3 2 2ay8 a 2 ) 

16 ir/t j-^sr - -^- + - j = 6^aX 2 . 



If X denote the force which must act upon the sphere in order to keep it at 
rest, we have 

XX = 6 ir/uaX 2 , 
or 

X = 6 iTixaX. 

These are the results given by Lamb in his treatise. In the general case the 
velocity u can be written in the form 



u 



-?[£{(« + 1> * + '£}- « (AM - * « WW] 
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with similar expressions for v and w. Applying the operator A 2 to this, we have, 
if we assume 

U = U x + U2 + .... + ttj, 



and 
also 



2 * * 

AX - % w 



*? u = ax I ' s i> 






The other terms vanish on applying this operator, as of course they should do. 

In the case where the axis of spin and the direction of the current velocity 
lie always in a plane normal to the surface of the solid, we have for the determi- 
nation of the pressure 



P 
p 



<*-(r+q) + k fff&JE**±* MjW . 



In the problem just discussed of the motion of the fluid all parallel to the axis 

of x, we have 

ug + vt] + w£ = 0, 

or the stream lines and vortex lines are at right angles to each other; this 

gives 

sin 8 = 1, 

and obviously in this case' the plane above mentioned is normal to the surface 
of the sphere. We have then 



P 
p 



G-(V+ 3 ) + lfff^<l*WM. 



Concerning Fwe have 

therefore 

V= 6 ir^aXx + const. 

In the general case, where there is no restriction as to the direction of motion of 
the fluid, the pressure must be determined by means of the equation 
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P 
P 



The computation would certainly be very difficult, if not impossible, from the 
complicated nature of the quantities involved. The first step, however, would 
be the determination of V; the angles 8 and $ can be found from the expres- 
sions for the velocities. If the motion is very slow and no external forces act 
on the fluid, we shall always have 



p 

p 



J (A 2 u . dx + A 2 v . dy + A 2 w . dz). 



Kirchhoff (vide Mathematische Physik, p. 377) has solved the problem of an 
ellipsoid of revolution rotating with constant velocity about its axis in a viscous 
fluid, both for the cases of an infinite extent of fluid and for a mass of fluid 
contained within a confocal ellipsoid. I do not see how to attack the general 
problem of the motion of any ellipsoid in a mass of viscous fluid ; but for 
the case of simple translation along one of the axes it is not difficult to find 
values for u, v, w which will satisfy all the prescribed conditions. Suppose an 
ellipsoid in the fluid with its axes coinciding with those of the co-ordinates, 

x l _i_ t + t = I 

If we assume first the case of no friction, and call V the potential of the 
ellipsoid at an external point, we have 

1 



V = nabc I — — d\l>, 

or 

V = const. - 2 77 (Ax 2 + By 2 + Cz 2 ), 

where cr is the greatest root of the equation 



a' + <r ' b* + <r ' c 2 + <r 



Now for the velocity potential <£ we have (American Journal of Mathematics, 
Vol. II. p. 260, et seq.), 

The quantities A, B, C are known to be 
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A = abc 



J (? 



d<jr 



u 

B = abcf — ^ , 

00 

C—dbcf *t ., 



in which 



JV= V(a 2 + ^) (6 2 + ^) (c 2 + V). 
The velocities u, v, w, in this case, will have for values 

- i i 1 ^ v 

U — — A-t- 2 »(2-^) 7fe»' 

i <? 2 r 

® ~" 2w{2- A) dxdy' 



w 



2tt(2—A) dxdz 

Reverting now for a moment to the case of the sphere, we had 

w = x _ 2Xl + |- W + g, 

u - % W + ^' 



in which 



.. \a? 1 

X 2 = X*r : 



a quantity proportional to the potential of the solid homogeneous sphere upon 

V 
an external point. The same remark, of course, holds concerning 2^(2 — A) ' 

Write now for the velocities in the case of the ellipsoid the following system of 
values similar to those obtained for the sphere : — 

, d , x . 1 d?V 

u = \-2V + Tx (xV) + 2wQi _ A) -&> 

d_ 1 <PV 

V ~ dy \ a ^ r ) "+" 2w (2 - ^) dxdy' 



to 



d_ 1 <Z 2 F 

<fo f^*) ~i~ 2tt (2 - .4) cfccfe' 
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These must satisfy the equation 

du . dv . dw 

di ~*~ % + Hz ~~ °» 



this gives 

-2£ + A»(**) = 0, 

or simply 

A 2 ¥ = 0. 

This is satisfied (Ferrer's Spherical Harmonies, p. 110) by assuming 

_ r df 

i/> and o- having the meaning already assigned them. For greater convenience 

we will write 

V= const. - (A 1 x 2 + B x f + Ciz 2 ), 
where 

A x = 2 it A, etc., 
also 

This (Ferrers, p. Ill) is the potential of a homogeneous ellipsoidal shell of 
determinate density at an external point. For the density we have (Kirchhoff, 
p. 179), 



which gives at once 



d<5? d& , 

dn { dn„ ' 



n U 4 ^ b 4 ^ c* 



that is, the density of the ellipsoidal shell at any point is proportional to the 
central perpendicular upon the tangent plane to the surface at that point. The 
values given now for u, v, w satisfy all required conditions, and we have in this 

case also 

Ahi . dx + A 2 « . dy + A 2 to . dz, 

an exact differential, viz. d. A 2 (x^). For greater generality, however, we may 
introduce two arbitrary constants, say a and fi ; then 

u = X - 2 a* + a Tx (*¥) + (4,-4) ^T> 



dy \°<W) "+" (4ir-^j) dbdy 

rf _J_ £ - d lZ. 

a dz \ x ^) + (4*- - A x ) dxdz 



v== a dy (^) + (4^-JO &*/' 



W 
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Now at the surface of the body we have 

u = v = w = 0, 
and also at the surface a- = ; therefore 

% = 2 irabc f v ^ + ^ ( / + ===. 

It will be convenient here to make a little digression and give the values of 
certain of our quantities as elliptic functions. Take \j/ 1} \jj 2 , xjj 3 as the variable 
parameters of a system of surfaces confocal to the given ellipsoid ; Sj, S 2 , 8 3 as 
the amplitudes of three elliptic integrals 

_ r <% a _ r &» a _ r <%* 

01 - J A (JcA) ' V* - J A (k&) > V*-J& (*& 



We have now 

, _ (a* + ft) (a 2 + ft) (a 2 + ft) 
* — (a? - b 2 ) (a 2 - c 2 ) ' 



, (5 2 + ft) (5 2 + ft) (6 2 + ft) 
2/ — (& 2 - c 2 ) (6 2 - a 2 ) ' 

(c 2 + ft) Q 2 + ft) (c 2 + ft) 
r — (c 2 - a 2 ) (c 2 - 6 2 ) 



Write now 



then make 



* 



.2 c* 



<*^Z*- t * 



£ = tan Si, 



where S t lies between and ■=. Similar transformations for ^ 2 and xp 3 give us 
finally 



/~2 -a- «£«• Ox en 2 en 6 3 

y — Va — <f ^ , 

on$ 1 sn0 2 dn 9 3 



z =Va*-<? m6i ■ 

In the above the modulus h is = \ a ~ „ • 

a? — c 2 

We have also for the quantities A x , B 1} C x the values 
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_ 2 aba */ l-dn»ft V 
■°» — (a 2 - c 2 )* J ^ dn 2 ^ J W l> 



_ 2«6c J sn 2 *, 







or 



^i "" ft 2 j K dk + ©(<?,) 



£ 2 = 



J (fl ^logiT 1 ©^-HT) ) 
A 2 j ft tfft + ^ © (flj + if ) j : 



where for brevity I have written 

2abc 

J- 775 „2\ 



(a 2 - c 2 )i 

The above transformations are given in full in an article, " On the Motion of 
an Ellipsoid in a Fluid," American Journal of Mathematics, Vol. II. We find, 
by the same transformations, 

n 7 C&h &irabC 

V = 2nabcJ£ i =^==0 1 , 
and also 

% = 2^abcJJ i =^ T ==K, 

K denoting the complete elliptic integral of the first kind. Equating now to 
zero the found values of u, v, w, we have at once, by means of the foregoing 
transformations, 

_ Q 2 - <?)i A. 

a ~ 2 abc ' tfj-l sn 2 1 dO 1 - 2 -k (a 2 - (?) K 
and 

fi = of a. 

The computation of the force necessary to keep the body in place, or the 
force which must be applied to overcome the resistance of the fluid if the body 
is moving through it, could be in this case, as in the former, computed by means 
of the dissipation-function ; but the process would be rather tedious and compli- 
cated. We can, however, without much difficulty, compute the pressure over 
the surface of the body supposed in motion. Observe that since ¥ is a surface 
potential corresponding to a surface density 



